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Optical measurement of the motion of a 940 kHz mechanical resonance of a silicon nitride nanostring res¬ 
onator is demonstrated with a read out noise imprecision reaching 37 dB below that of the resonator’s zero-point 
fluctuations. Via intensity modulation of the optical probe laser, radiation pressure feedback is used to cool and 
damp the mechanical mode from an initial room temperature occupancy of h\y = 6.5 x 10^ (r^ = 295 K) down 
to a phonon occupation of (n) = 66 ± 10, representing a mode temperature of Tm ~ 3 mK. The flve decades of 
cooling is enabled by the system’s large single-photon cooperativity (Ci = 4) and high quantum efficiency of 
optical motion detection (r|^ = 0.27). 


Cavity-optomechanical systems utilize multi-pass scatter¬ 
ing of light within a cavity to perform sensitive measure¬ 
ment of mechanical motion, with applications ranging from 
inertial microsensors HI 12 to transducers for interfacing dis¬ 
parate quantum systems dill. In this work, we integrate 
a silicon nitride nanostring mechanical resonator of motional 
mass meff = 90 picograms and frequency C0ni/27r = 940 kHz, 
with a photonic crystal optical nanocavity. The strength of 
the optomechanical coupling in this structure is characterized 
by a per photon measurement rate of the nanostring motion 
which is four times that of its intrinsic mechanical damping 
rate. Combined with an overall optical detection efficiency 
of r|^ = 0.27, this enables a measurement imprecision which 
reaches 37 dB below that of the zero-point fiuctuation noise 
of the bare mechanical resonator. Active cancellation of the 
mechanical thermal motion through feedback on the read-out 
laser’s intensity realizes cooling Eia from room temperature 
down to a phonon occupancy of (n) = 66 ± 10, corresponding 
to an effective mode temperature of Tm = 3 mK. This chip- 
scale microresonator, operating in a room temperature envi¬ 
ronment yet close to its quantum ground-state of motion, has 
a thermal-noise-limited force sensitivity of 125 aN/Hz^/^, a 
bandwidth of 190 kHz around resonance, and a linear dynamic 
range at one second integration time of greater than 60 dB. 

Resolved-sideband radiation pressure cooling has recently 
been demonstrated Em as an effective means to cool a tar¬ 
geted mechanical resonance of a structure close its quantum 
mechanical ground-state of motion. This technique, sharing 
similar physics to the resolved-sideband cooling of trapped 
ions requires spectral filtering of the upper (anti-Stokes) 
motional sideband from the lower (Stokes) motional side¬ 
band by a high-2 cavity in which the cavity linewidth (k) 
is narrower than the mechanical resonances frequency (0)m)- 
To date, experiments involving resolved-sideband cooling of 
mesoscopic mechanical objects to their quantum ground-state 
have relied on cryogenic pre-cooling using conventional re¬ 
frigeration means. Bath temperatures Tb < 100 mK are uti¬ 
lized for microwave devices m to enable high-2 supercon- 
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ducting cavities, whereas in the optical domain im more mod¬ 
est bath temperatures of ^ 10 K in a helium cryostat have been 
employed to, among other things, reduce intrinsic mechanical 
damping. 

An alternative method of radiation pressure cooling, one 
which is more amenable to lower frequency mechanical res¬ 
onators, relies on low noise optical read out of mechanical dis¬ 
placement combined with active feedback of the optical probe 
intensity Eiaiioi. Previous optomechanical feedback cool¬ 
ing experiments ifTTUT^ have demonstrated the ability to cool 
a wide range of mechanical resonators, from suspended large 
scale kilogram mass mirrors Ha to optically levitated micro¬ 
spheres EH [131. Although these experiments have realized 
substantial cooling, reaching the quantum mechanical ground- 
state of motion using active feedback cooling remains an illu¬ 
sive goal due to the stringent requirements on the measurement 
imprecision and the mechanical resonator 2-factor. Cooling to 
the ground-state requires a displacement measurement capa¬ 
ble of resolving motion at the level of the quantum zero-point 
fiuctuations of the mechanical resonator within its thermal de¬ 
coherence time, and with back-action close to the Heisenberg 
limit fT HTTljT^ . This regime has recently been approached in 
several cavity-optomechanical systems at liquid helium tem¬ 
peratures HIED, with feedback cooling of a MHz-frequency 
nanostring resonator being demonstrated in Ref. 021 down to 
an occupancy of 5 phonons. Here we set out to explore limits 
to feedback cooling of a photonic crystal nanostring structure 
in a room temperature environment suitable for a broad range 
of sensing applications. 

Cooling with simple derivative feedback can be understood 
from the following (classical) harmonic oscillator equation of 
motion, 

x+Yii + co^x= — - j, (1) 

OTeff V / 

where x{t) is the amplitude of motion of the mechanical res¬ 
onator, meff is the motional mass of the mechanical resonator, 
(Dm is the mechanical resonance frequency, yi the intrinsic me¬ 
chanical energy decay rate of the resonator, and TJh is the 
effective noise force of the thermal bath coupled to the me¬ 
chanical resonator. The final term on the right-hand side of 
eqn. O is the feedback forcing term, where g is the unit- 
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FIG. 1: a, False color SEM image of a device as used in the experiment, made from highly stressed, 415nm thick, stoichiometric silicon 
nitride released from a silicon wafer. The device consists of a photonic crystal cavity suspended on either side by nanoscale cross-section 
(115 ;umx 130 nm) tethers. The green and blue overlay regions indicate the photonic crystal cavity and the coupling waveguide, respectively. 
Due to the extreme aspect ratio of this device, only the ends of the tether region are shown, with the extent of the missing gaps on either side 
of the center of the device indicated, b, Zoomed-in SEM image of the photonic crystal section (green) and the adiabatically tapered on-chip 
coupling waveguide (blue), c, EEM-simulated electromagnetic energy density of the TE^ o optical resonance, with the outline of the Si 3 N 4 
beam shown in white, d, EEM simulation of the first-order in-plane differential mechanical resonance (IPj,o) of ^ dual nanobeam structure, 
indicating the displacement of the beams. Eor clarity, the tethers in the simulated structure are shorter than those of the actual device. In c 
the color scalebar indicates large (red) and small (blue) optical energy density, whereas in d the scalebar indicates large (red) and small (blue) 
displacement amplitude. 


less gain of the feedback loop, 0fb the phase of the feedback, 
y{t) = x(t) is the estimated resonator’s amplitude of 

motion from measurement, andxn(0 is the measurement noise 
(error). When Xn(^) is negligible and 0fb = 0, this term leads 
to viscous damping and cooling of the mechanical resonance, 
with the closed-loop mechanical Q-factor and phonon occu- 
pancy scaling as = 2m,i/(l +g) and (n) = nb/{l+g), re- 
spectively. Here, 2ni,i is the intrinsic mechanical quality fac¬ 
tor, n\) is the thermal bath occupancy at COm, and (n) is the 
resulting average phonon occupancy of the mechanical mode 
under feedback. For the room temperature measurements of 
this work, /Zb ~ = 6.5 x 10^, where 7b = 295 K 

is the bath temperature, is Boltzmann’s constant, and h is 
Planck’s constant over 2%. 

A more rigorous derivation of feedback cooling including 
quantum measurement noise O [T0]| shows that in the limit of 
large feedback gain (g ^ I) the cooled phonon occupancy is 
approximately given by, 

/ \ , ^ ^ f «b + l/2^ , f , nBA\ 

(”) + 2 ~ (- g - J + + — j ■ (2) 

Here zzimp is the open-loop measurement imprecision in units 
of photon number of the undamped oscillator, and zzba is 
the open-loop quantum back-action noise of the optical po¬ 
sition measurement. The first term on the right hand side of 
eqn. ^ represents the damped thermal noise from the bath, 
and is limited by the achievable feedback gain. The quan¬ 
tum fluctuations of the probe laser light manifest as radiation 
pressure shot noise ED, imposing a quantum-limited rela¬ 
tion between the imprecision and back-action noise sources, 
^BA,SN = 1 /(16r|tzzinip,SN), where r|t is the quantum efficiency 
of detection of light that enters into the optomechanical cav¬ 
ity and is scattered by the mechanical resonator. The last two 
terms on the right hand side of eqn. 0 thus represent a mea¬ 
surement limit to the attainable cooling. The optimal feed¬ 


back gain is gopt = \/(^b + ^ba) /^imp, which yields the min¬ 
imum attainable resonator occupancy assuming only quantum 
back-action, ((zz)niin + 1/2) = 2A^ni^^n\y A l/16r|t. Achieving 
{n) < 1 thus requires both zzimp < 1/2/Zb and r|t > 1/9. Con¬ 
sideration of the feedback bandwidth brings an additional con¬ 
straint m , requiring gni,i > ^b- In addition to these fundamen¬ 
tal constraints, technical limitations such as optical-absorption 
heating and thermo-refractive noise fT2\ may also play a role, 
adding excess back-action and imprecision, respectively. 

Shown in Fig. the optomechanical structure studied in 
this work consists of a “zipper” photonic crystal optical cav¬ 
ity 1^ supported by nanoscale tethers. The structure is fab¬ 
ricated using standard electron beam lithography and plasma 
etching techniques, and formed out of a 415 nm thick layer 
of stoichiometric silicon nitride (Si 3 N 4 ) deposited on a sili¬ 
con (Si) handle wafer. The zipper optical cavity (green shaded 
region of Fig. [^) consists of two micron-wide beams with lin¬ 
ear hole patterning, separated by a small gap of s = 150 nm, 
and attached to the bulk by 115 ^m long nanotethers of width 
w = 130 nm. The optical cavity design was simulated and opti¬ 
mized 1241 using the COMSOL flnite-element-method (FEM) 
mode solver. The optical cavity mode of interest is the fun¬ 
damental even mode (TE^ o) with electric field polarization 
predominantly in the plane of the Si 3 N 4 film. The central 
modification of the hole shape and location strongly confines 
the electromagnetic energy in the gap between the beams of 
the optical cavity (see Eig. [T]:), which results in a large shift 
in the TE^^o resonance frequency (cOc) with relative in-plane 
displacement of the beams. The mechanical mode of inter¬ 
est is the fundamental in-plane differential mode of the beams 
(IPd,o) depicted in Fig.[TJl. Eor the geometry considered here, 
the optical resonance wavelength is in the 1500 nm band 
(cOc ^190 THz), and the mechanical resonance frequency is 
0)ni/27r = 940 kHz. The simulated motional mass and zero- 
point motion amplitude of the IPj^o mechanical mode are 
^eff = 90 X 10“^^ kg and xzpf = 3 /^/(2meffOL)m) = 9.7 fm, 
respectively, where the generalized coordinate of mechanical 
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motion, x, is chosen to be the point of maximum in-plane dis¬ 
placement of the beams. 

To lowest order in the mechanical amplitude, the sensi¬ 
tivity of the optical resonance frequency to mechanical mo¬ 
tion is quantified by the linear dispersive coupling parame¬ 
ter goM = 9cOc/9v. For the optical TE^^o and the mechan¬ 
ical IPj 0 the optomechanical coupling is simulated to be 
= 41 GHz/nm, corresponding to a vacuum coupling 
rate of go = ^om-^zpf = 271 [358 kHz]. The interaction Hamil¬ 
tonian of the coupled cavity-optomechanical system is given 
by Hini = where he {rie) is the intra-cavity photon 

number operator (average). The corresponding average radia¬ 
tion pressure force applied by the intra-cavity optical field on 
the mechanical resonator is Tip = — (3^int/9v) = — %om^c- 
By modulating the laser intensity input to the optical cavity 
we can create a feedback cooling force as in eqn. Q. goM 
is therefore a critical parameter, determining both the optical 
measurement sensitivity to mechanical motion and the strength 
of the radiation pressure force that can be applied per photon. 
In principle, small goM can be overcome with larger optical 
power. In practice, parasitic effects such as thermo-refractive 
noise 12^ and optical-absorption heating then tend to limit the 
achievable measurement imprecision and back-action ifT^ . In 
this regard, a figure of merit is the single-photon cooperativity 
(Cl), which physically represents the ratio of the displacement 
measurement rate per photon to the decoherence rate per (ther¬ 
mal bath) phonon of the mechanical resonator. For the zipper 
cavity studied here, Ci = 4gQ/Kyi = 4. 

The experimental setup used to characterize the optical 
and mechanical properties of the zipper cavity is shown in 
Fig. [^a). The device is mounted inside a vacuum chamber 
reaching a pressures of 2.5 x 10“^ mbar, sufficient to eliminate 
the effects of gas-damping of the mechanics. A tunable exter¬ 
nal cavity semiconductor diode laser (New Focus Velocity se¬ 
ries) is used to provide both the signal beam on-resonance with 
the zipper cavity and the local oscillator beam (LO) for homo¬ 
dyne detection. Laser light is efficiently coupled into and out 
of the zipper cavity using an optical fiber taper 1^ in combi¬ 
nation with an on-chip tapered waveguide (blue shaded region 
of Fig. [^. Tapering the width of the on-chip waveguide allows 
for adiabatic mode-conversion between the waveguide and the 
tapered fiber placed upon it as in Ref. ll^ . The reflected opti¬ 
cal signal from the zipper cavity is separated using an optical 
circulator and sent to a balanced homodyne detector (BHD). A 
low-pass filtered (LPF; bandwidth < 200 kHz) version of the 
BHD signal is sent to a control circuit (PID) which drives a 
fiber stretcher (FS) to lock the relative phase (0/^) between op¬ 
tical LO and reflected signal beam, and sets the phase quadra¬ 
ture of the homodyne detected signal. A band-pass filtered 
(BPF; bandwidth 0.2-1.9 MHz) version of the BHD signal is 
sent to an electronic spectrum analyzer (ESA) to measure the 
mechanical noise spectrum. In the case of optical feedback, 
the band-pass filtered signal is also sent to an analog differen¬ 
tiator circuit (Toptica mFALC), whose output is sent through 
a variable phase shifter, A(|), and finally onto an electro-optic 
intensity modulator (IM) which closes the feedback loop and 
modulates the signal beam intensity. 

In Fig.[^ we plot the reflected optical power normalized by 



feedback circuit 



FIG. 2: Device Characterization and Experimental Setup, a, 

Schematic of the optical and electrical set-up used to characterize 
and feedback-cool the mechanical resonator. Additional acronym: 
fiber polarization controller (FPC). b. Plot of the laser reflection 
spectrum when scanned across the optical mode used in the exper¬ 
iment (blue) and fit (red). Dashed grey line indicates on-resonance 
laser position during measurements. The measured loaded optical 
g-factor is Qq = 2.5 x 10^, with a waveguide loading to total cavity 
loss rate ratio of Ke/K = 0.63. c, Power spectral density (PSD) of 
the mechanical resonator’s thermal noise near 940 kHz, transduced 
using the setup in (a) with the laser on-resonance, A = 0, and at low 
power ^ 1 • d» Autocorrelation of the slowly varying amplitude of 
the transduced mechanical thermal noise, (x) = {\x{t) | |x(7 -hx) |). 
x{t) = x(7) sin(cOmt + 9(0)) defines the slowly varying amplitude, 
x{t). An exponential fit (red curve) to the measured data (green cir¬ 
cles) yields an intrinsic mechanical damping rate yi/27i = 1.76 Hz 
and quality factor j = 5.3 x 10^. 


the input power at the device as the laser frequency is scanned 
across the TE^ o optical resonance of the zipper cavity. The 
background level of this normalized plot, Tjr = 0.43, deter¬ 
mines the overall detection efficiency of light emitted from 
the cavity, including the single-pass coupling between the chip 
waveguide and optical fiber taper (r|c = 0.72), optical loss be¬ 
tween the fiber taper and BHD photodetectors (r|f = 0.69), and 
the quantum efficiency of the photodetectors (tiqe = 0.88). 
The measured linewidth of the optical cavity resonance is 
k/2tz = 11 GHz, corresponding to a loaded optical g-factor 
of gc = 2.5 X 10^. A low optical g was chosen to increase 
the optical power handling ability of the zipper cavity. From 
the depth and Fano-like shape of the resonance dip, the cou¬ 
pling rate of the cavity to the on-chip waveguide is estimated 
to be Ke/27r = 49 GHz, yielding a slightly overcoupled system 
with Ke/K = 0.63. The total efficiency of motion detection is 
related to the detection of laser light that enters the cavity and 
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interacts with the mechanical resonator, which for this device 
and set-up is r|t = r|r(Ke/K) = 0.27. 

Measurement of the mechanical motion, x{t), is performed 
by setting the laser frequency (coi) to the resonance frequency 
((Dc) of the optical mode (A = coi — (Dc = 0), and monitoring the 
phase quadrature = 7r/2) on the balanced homodyne detec¬ 
tor |[2Q| . The measured single-sided noise power spectral den¬ 
sity (NPSD) is plotted without feedback (open loop) in Fig.|^, 
showing the optically transduced thermal Brownian motion of 
the IPj^o mechanical resonance at 940 kHz. Here, the NPSD is 
calibrated using the measured goM of the IPj^o mechanical res¬ 
onance (along with optical input power, optical cavity param¬ 
eters, and optical detection efficiency), and plotted as Syy{(Si) 
in units of m/\/Hz. The optomechanical coupling is inferred 
from the optical spring shift of the mechanical resonance fre¬ 
quency 1^ . and is measured to be goM/^TT = 36 GHz/nm, in 
close agreement with the simulated coupling coefficient for the 
measured device geometry. Autocorrelation of the magnitude 
of the slowly varying amplitude of the resonator thermal mo¬ 
tion, ^|jc|(x) = (|v(r-|-x)||v(r)|), is shown in Fig. 1^. A fit to 
the exponential decay of (x) (red curve) yields an intrinsic 
mechanical damping factor of yi/27r = 1.76 Hz, corresponding 
to an intrinsic mechanical Q-factor of Qm,i = 5.3 X 10^ 

Figure shows an open-loop, wideband spectrum of the 
measured mechanical NPSD (A = 0, 0/^ = 7r/2) for an in¬ 
tracavity photon number of Hq = 0.17 (blue curve) and for 

= 0 (grey curve). Measurement of the noise level with 
the signal arm blocked corresponds very closely to the sig¬ 
nal vacuum-noise level (electronic detector noise is 12.9 dB 
below the measured noise level). Here we have plotted the 
NPSD in units of phonon quanta, Syy{(o) /(•x|ppyi/4). In these 
units the mechanical mode occupancy {{n)) and the phonon 
imprecision level (^imp) can be simply read off from the peak 
height of the 940 kHz resonance and the nearby background 
level, respectively HSl. A plot of the measured ^imp versus 
Uc is shown in Fig. as grey circles. The expected impre¬ 
cision due to vacuum noise of the signal beam is ^en by 
^imp,SN = i^yi/(64^c^ofit)’ which is plotted in Fig. ^ with 
no free parameters for the measured r\t = 0.27 (solid cyan 
curve). Also plotted are the theoretical quantum back-action 
due to the shot noise of the signal beam (^ba,sn; solid red 
curve) and the quantum-limited total added measurement noise 
(solid green curve). The minimum total added measurement 
noise assuming ^ba = ^ba,sn occurs at a signal power corre¬ 
sponding to He = 0.12, and represents the standard quantum 
limit (SQL) for our measurement set-up, ^sql = 1/(2yfifi) = 
0.96 quanta i27l . The imprecision is vacuum noise limited 
for all but the highest powers (Uc > 500), reaching a value 
34 dB below the SQL imprecision for an ideal detector (=1/4 
quanta ll28l ). For all the measurement powers shown, the back- 
action noise is bounded (^ba < 4 x 10^) by the standard devi¬ 
ation (±3%) in our measurement of the large thermal noise in 
the mechanical resonator. 

Feedback cooling of the main 940 kHz mechanical reso¬ 
nance is performed as follows. As in the open loop measure¬ 
ments, the laser is tuned on resonance with the optical cavity 
(A = 0) and the LO phase is locked to 0/^ = 7r/2 which max¬ 
imizes the measured BHD signal due to mechanical motion. 




FIG. 3: Measurement noise, a, Measured wideband NPSD plotted in 
units of phonon quanta of the dominant mechanical mode at 940 kHz 
for He = 0.17 (blue) and with the signal beam blocked (ric = 0; grey 
curve). The insets of (a) are simulated mechanical mode displace¬ 
ment profiles for the two mechanical resonances visible in the spec¬ 
trum. The mode at 740 kHz is the poorly transduced in-plane com¬ 
mon motion of the zipper beams. Frequencies below 200 kHz are 
attenuated by a high-pass filter and the small slope of the background 
noise level is a result of the frequency-dependent gain of the balanced 
photodetectors, b, Plot of the measured imprecision noise in units of 
phonon quanta (^imp; grey filled circles). Also plotted are the theo¬ 
retical vacuum-noise-limited imprecision (wimp,SN^ solid cyan curve), 
theoretical quantum back-action noise (n^A,SN'^ solid red curve), and 
theoretical quantum-limited measurement noise (wimp,SN + ^BA,SN 5 
solid green curve), assuming the measured r|r = 0.27 and with no 
additional fit parameters. Dashed curve is the measurement impreci¬ 
sion for an ideal continuous position measurement with = 1. 


The BHD signal (y{t)) is fed to an electronic feedback circuit 
which modulates the intensity of the probe laser incident on 
the cavity approximately in proportion to —y. The mechani¬ 
cal resonator responds to the intensity modulations of the inci¬ 
dent probe laser, imprinting its motion in the orthogonal phase 
quadrature of the reflected probe light for A = 0 1^ . As de¬ 
sired then, for this configuration (A = 0, 0/^ = :7r/2) it is only the 
response of the mechanical resonator to the feedback which is 
recorded in the BHD signal. The electrical gain of the feed¬ 
back circuit is held fixed at a value found to yield maximum 
cooling at the highest value of ric (=734), and the laser probe 
power is increased from low to high, increasing the total loop 
gain, and thereby increasing the observed cooling and damp¬ 
ing. 

Figure El shows the measured NPSD around the main 
940 kHz mechanical resonance with optical feedback applied 
as per the above prescription. At each measured optical power, 
the mechanical spectra are fit with a Lorentzian curve (black 
solid lines) from which an area and linewidth of the spectrum 
are determined. The phonon occupancy of the mechanical res¬ 
onance, plotted as blue circles in Fig. |^, is proportional to 
the transduced area under the mechanical spectrum normalized 
by He, whereas the damped mechanical Q-factor is determined 
from the linewidth. Absolute calibration of the phonon occu¬ 
pancy for each of the optical powers is found by comparing to 
the measured area under the mechanical spectrum at the lowest 
power point (Uq = 0.36) with the feedback off. At this power, 
and in open loop, dynamic back-action effects are negligible 
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FIG. 4: Room-temperature laser cooling, a, Plots of the transduced 
mechanical spectra with the optical feedback engaged for increasing 
Uc, showing the damping and cooling of the dominant mechanical 
mode. Color scale indicates photon number of each mechanical spec¬ 
trum. The black curves are the fits to the measured spectra used to 
extract the damped mechanical g-factor (Qm) and phonon occupa¬ 
tion ({n)). Spectral peaks at 710 kHz and 910 kHz are weakly trans¬ 
duced mechanical modes, while that at 1 MHz is extraneous elec¬ 
tronic noise. Inset: room-temperature mechanical spectra taken with 
the feedback off for calibration, b, Plot of the inferred (n) (blue cir¬ 
cles; left axis) and gm (green solid circles; right axis) extracted from 
the fits to spectra shown in (a). The first data point (square), at low¬ 
est He, is taken from the data in the inset of (a) with the feedback off 
and is used as a room-temperature calibration. The y-axes scales are 
normalized such that ideal cold-damping would result in the blue and 
green data points lying directly on top of each other. Solid black curve 
is the theoretical prediction of cooling with no fit parameters and the 
dashed line is the minimum possible cooling for our system if the cir¬ 
cuit gain is optimized at each ndseQ App. 0. The uncertainty in the 
inferred (w) of ±15% is smaller than the data points and is dominated 
by the 95% confidence interval of the fits to each mechanical spec¬ 
tra. Inset shows the measured phonon occupation versus mechanical 
quality factor during the cooling run (blue points). A linear fit to the 
data is shown as a red curve with R squared value of 0.9988. 


and the mechanical resonance is at the room temperature ther¬ 
mal occupancy (Tb = 295 K, (n) = fib = 6.5 x 10^). Of cru¬ 
cial importance to the interpretation of the data presented in 
Fig. 1^ is the fact that the change in the mechanical quality 
factor (green dots) follows the change in the measured occu¬ 
pation (see inset), which as discussed earlier, is a hallmark of 
feedback cooling. 

The lowest phonon occupation achieved here is (n) = 66 ± 
10. At this cooling point the linewidth of the mechanical res¬ 
onance at C0ni/27r = 940 kHz has been broadened to 190 kHz. 
Further cooling of the mechanical resonance is limited primar¬ 
ily by the combination of two sources. An excess time delay 
of l^s in our feedback loop modifies the broadband phase re¬ 
sponse of the system away from the ideal value of 0fb = 0, 
leading to amplification, rather than damping, of mechanical 


motion for frequencies outside an approximate 250 kHz band¬ 
width around the mechanical resonance (see App. [^. Addi¬ 
tionally, at the highest optical cooling power of ric = 734 we 
are nearing an optical-absorption-induced, thermo-optic bista¬ 
bility of the optical cavity response. Calibration of the absorp¬ 
tion heating via the thermo-optic tuning of the optical cavity 
resonance indicates only a AT 10 K increase of the local bath 
temperature for = 734, representing a parasitic back-action 
noise of ^BA,abs ~ 300^c- Although ^BA,abs ^ ^b at this optical 
power, which doesn’t impact the current level of cooling, as 
one approaches mode occupancies of {n) < 1 the back-action 
noise becomes relevant. 

Cooling to the quantum ground-state from room temper¬ 
ature remains an achievable goal, but requires improvement 
in several key device properties. As discussed in Ref. ISl, 
one cannot increase the feedback bandwidth without limit as 
eventually the amount of imprecision noise (white shot noise 
in the ideal case) fed back onto the mechanical resonator 
is enough to heat it out of the ground-state. As shown in 
App. E consideration of the feedback bandwidth in turn con¬ 
strains the intrinsic mechanical g-factor to a value greater 
than the thermal bath occupancy (Qm ^ 3fib). Quality fac¬ 
tors approaching this limit for MHz-frequency resonators in 
thin film Si 3 N 4 have been achieved by modifying the structure 
to minimize losses at the clamp points (301 [JU, utilizing ni¬ 
tride with fewer bulk defects and higher stress (321, and mod¬ 
ifying the post-etching surface properties of the nitride Ea. 
Even with an increase in the mechanical g-factor, however, 
reaching the ground-state with the optomechanical coupling 
strength of the devices in this work would still require a pro¬ 
hibitively large intra-cavity photon number of ~ 5 x 10"^. 
By increasing the optomechanical coupling and optical g- 
factor to levels previously demonstrated in similar devices 
(^OM/27r = 200 GHz/nm, g^ = 6 x 10"^ (34l), this can be low¬ 
ered to a value as small as Hq ~ 60. There is also the issue of 
the parasitic back-action. In the current devices we estimate an 
imprecision-back-action product of ^imp^BA ~ 450(^inip^BA)H 
at He = 734, where (^imp^BA)H = 1/16 is the Heisenberg quan¬ 
tum limit. With the above mentioned device improvements, 
^BA,SN ^ ^BA,abs 5 und the imprecision-back-action product 
would approach l/r|t ~ 3.7 of the Heisenberg limit, suitable 
for cooling below unit occupancy. 

Viewed from the perspective of a continuous position mea¬ 
surement, attaining ground-state cooling requires (approxi¬ 
mately) that the rate at which measurement information is 
gained about the mechanical motion (Fmeas = 4r|^^c^o/^) 
greater than the rate at which the mechanical resonator is dis¬ 
turbed by coupling to its thermal bath (Eth = yi(fib ± 1)), and 
approach the back-action decoherence rate (F^ = Yi^ba) (FTI 
QiiEa. More generally then, the photonic crystal optome¬ 
chanical devices studied here could enable quantum measure¬ 
ment and control protocols (36ti4^ for the preparation of me¬ 
chanical objects in highly non-classical quantum states of mo¬ 
tion. With the device improvements mentioned above, these 
protocols could be implemented without additional cryogenic 
cooling, and in a room temperature environment where they 
may be employed for a variety of precision sensing appli¬ 
cations. In a similar vein, feedback control is commonly 
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employed in MEMS sensors of forces and fields B3l 1441 to 
change the frequency, bandwidth, and dynamic range of the 
sensor. In this work the bandwidth of the resonator’s response 
is increased from 1.7 Hz to 190 kHz and the dynamic range is 
increased by 50 dB, all while preserving the undamped ther¬ 
mal noise force sensitivity of 125 aN/Hz^/^. When applied to 
the field of atomic force microscopy, for example, such an op¬ 
tomechanical sensor isi might be used to improve imaging 
resolution by reducing thermal motion of the sensor tip 1461 . 
or in the case of measurements of molecular motion, to resolve 
dynamics at microsecond time scales El. 
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FIG. 5: Normalized optical reflection spectrum of the device studied 
in the main text, measured by scanning the probe laser wavelength 
from X = 1540-1570 nm. The signal level is the ratio of the measured 
reflected optical laser power to the optical power incident at the opti¬ 
cal cavity. The blue curve is the measured laser wavelength scan and 
the red curve is a fit used to determine the optical coupling efficiency 
between the collection flber and the on-chip waveguide, as well as the 
optical decay rate (k) and the waveguide-to-cavity coupling rate (Ke) 
of each optical cavity mode. The images below the normalized opti¬ 
cal reflection scan are the numerical simulations of the energy density 
of the identifled optical modes. The simulation images are labeled as 
TEg('^) where the subscript e{o) indicates the spatial mode symme¬ 
try about the Ox mirror plane and the subscript indicates the order of 
the mode along the waveguide axis. 
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O. Painter, Optics Express 17, 20078 (2009). 
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A. Device characterization 
1. Optical mode characterization 

The photonic crystal cavity used in this work exhibits a se¬ 
ries of optical resonances of different symmetries and orders. 
A detailed description of a photonic crystal with a similar 
structure can be found in Ref. |[34l . The main symmetry of 
interest (Cx) is about a plane which runs parallel to the length 
of the beams and is situated in the center of the gap between 
them. The optical field spatial mode pattern can be either even 
(e) or odd (o) about this symmetry plane, where e modes have 
their fields concentrated in the gap and the field of o modes 
experience a node at the center of the gap. Due to their higher 
field concentrations in the gap, the resonance frequencies of e 
modes are more sensitive to the in-plane differential motion of 
the two beams, x, than are the o modes. As a result the e modes 
tend to have larger optomechanical coupling to the in-plane 
differential motion of the nanobeams. By measuring how well 
an optical mode transduces differential motion of the beams 
we can determine its symmetry. Accordingly, the resonances 



PIG. 6: Measurement of the mechanical spring shift for a range of 
laser-cavity detunings (blue circles). The theoretical frequency shift 
for an optomechanical coupling coefficient of goM/^Tt = 36 GHz/nm 
(go/2n = 358 kHz) is shown as a red curve. 


at X = 1551 nm and 1568 nm are labeled as odd symmetry, 
whereas the optical mode of interest at = 1561 nm is of even 
symmetry. In addition the optical modes are labeled by an in¬ 
teger subscript 7 = 0 , 1 , 2 ,... related to the number of nodes in 
the spatial mode pattern of the optical mode along length of 
the beams. The order of the optical modes are determined by 
comparing the measured optical frequency to that found from 
COMSOL Multiphysics FEM simulations ||48]| . 

Figure shows the laser reflection spectrum of the device 
“zipper” cavity device studied in the main text. There are three 
optical resonances measured in the 1550 nm telecom wave¬ 
length band, which we identify and label as in the simulated 
plots shown below the measured reflection spectrum in Fig.|^ 
The resonance at X= 1561 nm is the one used for feedback 
cooling in the main text. Considering the Fano-like lineshape 
of the resonances, we fit the wide optical scan (blue curve) to a 
model consisting of three overlapping optical resonances and 
find good agreement (red curve). The fit to the TE^ 0 mode 
of interest at 1561 nm yields (Oc/27r = 192 THz, Qc = 2500, 
K/2;r = 77 GHz, and Kg /k = 0.63, where cOc is the angular fre¬ 
quency of the optical cavity, Qc is the quality factor, K is the 
total angular loss rate, and Kg is the angular coupling rate to 
the waveguide. 

2. Calibration of the optomechanical coupling strength 

In order to calibrate the optomechanical coupling coeffi¬ 
cient, goM = d((ic/dx, we measure the optically-induced spring 
shift of the IPj^o mechanical resonance versus laser-cavity de¬ 
tuning A. In the sideband unresolved regime in which (Om <C K 
the mechanical frequency shift is given by 1 ^ . 


COm 




/ \ A 

V Weff< J A2+(|)^ 




%OM^c A 


A 

a^ + (!A 


(Al) 


where 0 )ni,o is the bare mechanical frequency, ric is the number 
of intracavity photons, meff is the effective motional mass, K is 
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FIG. 7: Measurement of the amplitude autocorrelation, {\x{t)\\x{t-\- 
t)|), of the thermal motion of the IPj o mechanical resonance. The 
blue data is the measured autocorrelation and the red line is a fit to the 
exponential decay, yielding a mechanical g-factor of gm = 5.3x10^. 

the optical cavity (energy) decay rate. In Fig. we show the 
measurement of the spring shift (blue data) and the fit curve 
(red line), from which we find the vacuum coupling rate to be 
go = goM-^zPF = 2;r[358 kHz]. The zero-point fiuctuation am¬ 
plitude vzpF = (^/2meffC0m)^^^ = 9.7 fm of the IPj^o mechan¬ 
ical mode is determined from the mechanical frequency and 
the simulated motional mass, m^ff = 90 x 10“^^ kg 14^ . This 
yields a linear dispersive coupling coefficient of goM/^TT = 
36 GHz/nm from the vacuum coupling rate, in good agreement 
with the simulated value from COMSOL ||48]| of goM/^TT = 
41 GHz/nm. 

3. Measurement of mechanical quality factor 

Here we describe details of the measurement of the intrin¬ 
sic mechanical quality factor of the fundamental in-plane dif¬ 
ferential mechanical resonance (IPj,o) of the “zipper” cavity 
structure. The approximately 1 Hz linewidths of this mechan¬ 
ical resonance combined with the minimum resolution band¬ 
width of our spectrum analyzer (also 1 Hz) makes measuring 
the mechanical quality factor from the spectral density diffi¬ 
cult. However, the amplitude of a thermally-driven resonator 
will be correlated with itself over a timescale dictated by the 
mechanical damping rate corresponding to the coupling of the 
resonator to the thermal bath. Formally it can be shown that 
the autocorrelation of the amplitude, {\x{t) \ |x(r H-x) |), will de¬ 
cay as where Xm = and (Dm is the intrinsic 

mechanical angular frequency lISOl . The measurements of the 
thermal amplitude of motion (v(r)) are made using the exper¬ 
imental setup described in the main text with the laser on- 
resonance with the optical cavity and the feedback off. The 
slowly varying envelope of the mechanical resonator motion, 
|x(r)|, is obtained from the magnitude channel of a lock-in 


amplifier tuned to the mechanical resonance frequency. We 
use a bandwidth of ~ 100 Hz, much larger than the mechan¬ 
ical linewidth, which ensures that any frequency jitter in the 
mechanical resonance does not affect the measurement. Fig¬ 
ure |7] displays the autocorrelation of the measured signal. A 
exponential curve fit to the decay of the autocorrelation signal 
yields a coherence time of Xm = 90 ms, which for (Dm/27r = 
940 kHz gives gm,i = 5.3 x 10^. For lower-g structures it was 
confirmed that this technique agrees with a direct measurement 
of the linewidth as measured on a spectrum analyzer. Fur¬ 
thermore, to be sure that we are not doing significant dynamic 
back-action, and thus that we have obtained the intrinsic qual¬ 
ity factor, we performed the measurement at the low optical 
power where there was no power dependence to the measured 

^m- 

B. Optical feedback cooling: theory 

In this section we present (without derivation) the results 
from Ref. 0 for the variances of the position and momentum 
of a harmonic oscillator in a derivative feedback loop with an 
on-resonant quantum noise limited laser. First, though, we de¬ 
fine the variables used in Ref. i): 


go = gOM-XZPF, 

(Bl) 

G = go\/^, 

(B2) 

K/i 

«imp,SN 

(B3) 

4G^ 

«BA,SN = —— , 

KYi 

(B4) 

_ 4gcdG0)m 

8 : 

(B5) 


where G is the parametrically-enhanced coupling rate between 
the optics and the mechanics, n^ is the average intracavity pho¬ 
ton number of the probe laser, ^imp,SN is the shot-noise-limited 
imprecision in units of phonon number, ^ba,sn is the quantum 
back-action of the shot noise in units of phonon number, r|t is 
the total quantum detection efficiency |[27l, g is a normalized 
unitless feedback strength, equivalent to that used in the main 
text, and gcd is a unitless gain term accounting for the feedback 
circuit response. 

The feedback response is taken to be a standard deriva¬ 
tive high-pass filter with cut-off frequency (Ofb, which in the 
Fourier domain is given by IT((d) = —/cogcd/(1 — 

With this assumed feedback response function the variances 
of the two quadratures of the mechanical mode are given by 
(under certain assumptions 0 valid in this work), 
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g ^imp + ( + ^BA 1(1“^- Y 
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(B7) 


To determine the fundamental cooling limits from these 
equations we take the limit of large feedback bandwidth 
(o)ni/o)fb)^ ^ 1 and large feedback strength g ^ 1, which al¬ 
lows us to drop some terms, and write a simpler formula for 
the position fluctuations, 


/s 2\ _ , «ba , nb + 1/2 

\Oq ) — g^imp H“ "k 

g g 


(B8) 


Taking the same limits for the momentum variance we And, 


(5/) = 


g^imp ■ 


^BA 
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^b +1 /2 

g 


H-Yi 


g ^imp0)fb 1 /_ X 

-2-^-(^b +^Ba) 

ff>m ®fb 


(B9) 


Here we note that the first bracketed term on the RHS of the 
formula for the momentum variance looks the same as that for 


the position variance in eqn. (B8). For now we will assume that 
the second bracketed term on the RHS is small and ignore it. 
We will revisit this assumption at a later point and determine 
when this assumption is valid. To be explicit, our working 
assumption is stated below: 


Yi 


g ^imptt)fb 1 /_ X 

—— + — (^b + ^ba) 

0)m t^>fb 


< 1 . 


(BIO) 
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' g^imp ■ 


^BA 
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(B14) 


The optimal g that minimizes (n) in eqn. (B14) is gopt = 

[(flb + ^ba) ^imp] which upon re-substitution yields a min¬ 

imum phonon occupancy of, 

(^min) 2 ~ ^imp (^b ^Ba) (B15) 

’ (B16) 


'2\ ^imp^b ■ 


16%’ 


where in the final expression we have used the quantum- 
limited relation, ^ba,sn = 1 /(16T|t^imp,SN)- From this final re¬ 
sult we can establish the two requirements mentioned in main 
text for reaching (n) <1, 


1 

At >9, 


^imp 


J-U-L 

16 ^b V At 


(B17) 

(B18) 


This leaves us with the simplified (and approximate) form of 
the momentum variance under optical feedback cooling. 


Of course other technical limitations may come into play, such 
as excess back-action or instabilities which limit the optical 
power and thus feedback gain. Furthermore, we must be in a 
parameter regime that satisfies the assumption we made at the 


outset in eqn. (BIO). Substituting gopt into eqn. (BIO) yields 
the following relation. 


/S 2\ , «BA , «b + l/2 

{op ) =gnimpH-^-• 

g g 


(BID 


2m,i > (^b “1“ ^Ba) ( Y 
\ tOm / 


(B19) 


In order to relate these variances to a phonon occupation 
number, we equate the total energy of the oscillator to the sum 
of its variances. 


E = [(S^^) + (5/)] = ^(n) + T/j (B12) 

which yields for the average mode occupation, 

{n) + l = \[{5q^) + {5p^)]. (B13) 

Substituting the simplified formulas for the variances 
(eqn. ( |B8| ) and eqn. ( |B11| )) gives eqn. (2) of the main text. 


where we have again used ((Om/cOfb)^ ^ 1. This relation can 
be further simplified by saturating the inequality in eqn. ( |B18 ) 
and relating ^imp to ^ba through their quantum-limited rela¬ 
tion. 


2m,i>(nb+9^Y-i)(^)' 

If we further assume ideal detection (r|t = 1) and take the feed¬ 
back bandwidth to be 0L)fb,opt ~ 3(Dm, a value numerically found 
to be the optimal feedback bandwidth in most cases O, this 
yields the requirement stated in the main text: 2m,i ^ 3/Zb in 
order to reach a cooled phonon occupation of (n) < 1 using 
optical measurement and radiation pressure feedback. 































10 



FIG. 8: a, Color density plot of {n) for ideal radiation-pressure derivative feedback cooling using the device parameters of the main text, versus 
normalized gain and photon number. The black line is a contour at {n) = 66, the minimum phonon number achieved in the main text. The 
white dashed line denotes Hq = 734, which is the maximum number of photons that can stably occupy the cavity. The minimum occupation 
is (^min) = 64 for Uc = 734. The global minimum at the ideal power and gain is (^min) = 8. b. Same as (a) except with perfect detection 
efficiency r|t = 1- The minimum occupation is (^min) = 38 for ric = 734. The global minimum at the ideal power and gain is (^min) = 6 . c. 
Same as (a) except with gni,i = 1-5 x 10^. (^min) = 9.4 for ric = 734 and (^min) = 1-8 for the ideal power and gain, d. Same as (c) except with 
Tjt = 1. (wmin) — 5 for Uq = 734 and (^min) = 0-95 for the ideal power. 


We can also use the full formulas, eqn. ( |B6| ) and eqn. ( [B7| ), 
to predict the maximum cooling for a given set of param¬ 
eters. Figure shows the theoretical cooling plot for the 
parameters of the device studied in this work. Figure 
shows the same cooling plot but for unit detection efficiency 
(r|t = 1). From these plots we see that for a maximum intra¬ 
cavity photon number, = 734, the best achievable for the 
device parameters in work and assuming ideal derivative feed¬ 
back cooling is (^min) = 64, very close to the measured mini¬ 
mum mode occupancy of {n) = 66 ± 10. With the limited de¬ 
tection efficiency of r|t = 027, the minimum achievable occu¬ 
pation is (^min) = 8. For perfect detection efficiency (^min) = 6 
could be reached. Both of these estimates have neglected ex¬ 
cess back-action. Furthermore, achieving this level of cooling 
requires using very high intracavity photon numbers, which 
would cause thermal instabilities in our system. As discussed 
in the main text, several device improvements would be re¬ 
quired to achieve this level of cooling. Figures. and d show 
cooling plots assuming an improved mechanical Q-factor of 
= 1.5 X 10’. 


C. Detection and feedback calibration 
1. Determination and locking of homodyne phase 

The relative phase angle between the local oscillator (LO) 
and the reflected cavity signal (0/^) determines which optical 
quadrature of the reflected signal one measures at the out¬ 
put of the BHD. One can easily show that for an optical 
probe on-resonance with the optical cavity (A = coi — COc = 0) 
the position of the mechanical resonator is imprinted on the 
phase quadrature of the reflected light corresponding to 0/^ = 
7t/2 1^ . This phase angle is set by the relative path lengths 
of the LO and signal arms, and slow thermal drifts or acous¬ 
tic noise can cause changes in the relative path length, thus 
changing the phase angle and the measured optical quadra¬ 
ture. Note that to avoid complications of laser phase noise, 
the LO and signal beams are derived from the same laser and 
the two optical paths are roughly matched in length. In order 
to stabilize % the low-frequency component of the BHD sig- 
nal (/ < 200 kHz) is fed to a digital feedback system (Toptica 
DigiLock) that has a high-voltage output panel (Vout < 150 V) 
which drives a voltage-controlled flber stretcher (Optiphase 
PZl). This flber stretcher (FS) consists of a long distance of 
optical flber 10 m) wrapped around a bulk piezo element. 
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FIG. 9: Plot of a balanced homodyne signal (blue) as the fiber 
stretcher voltage is swept (green), with the laser on-resonance with 
the optical cavity. 

Shown in Fig.|^is the homodyne output while linearly driv¬ 
ing the FS and for optical probing at A = 0. This voltage sweep 
(green) modulates the detected phase quadrature of the re¬ 
flected light, which we detect as a DC signal with sinusoidally 
modulated voltage output proportional to cos0/^. The thermal 
Brownian motion of the strongly coupled mechanical mode at 
/m = 940 kHz is evident in the blue curve as increased noise on 
the signal, which decreases rapidly at the maximum and min¬ 
imum points of the curve corresponding to the read out of the 
intensity quadrature of the reflected optical signal. Read-out 
of the phase quadrature of the reflected light (0/^ = %I2) oc¬ 
curs at the average voltage between the peak and trough of this 
curve (grey dashed line), which is the lock-point desired in our 
measurements. In the cooling curve shown in the main text, 
for each optical probe power the laser-cavity detuning and the 
homodyne lock-point are adjusted before each each mechani¬ 
cal spectrum is measured, and then checked after the measure¬ 
ment is completed and before moving to the next optical power 
point. 

2. System response and delay times 

Here we report the system response and delay times of our 
feedback circuit used to apply radiation pressure feedback. In 
Fig. we present the measured amplitude, phase response, 
and time delay of the analog circuit, mFALCllO, used as 
part of the feedback loop to perform the derivative feedback. 
The cyan dashed line in Fig. [T^ indicates the amplitude re¬ 
sponse for purely derivative feedback, which is proportional 
to CO. Over the frequency range used for this experiment (0.8 
- 1.2 MHz), the amplitude response deviates from the ideal 
slope by 0.5 dB and the phase changes hy ^ 5^. The time 
delay from the mFALCllO circuit, shown in Fig.p^, is neg¬ 
ligible when compared to the mechanical period. 

In Fig. we present the amplitude and phase response of 
the entire feedback loop, including the optical portion. It is im¬ 
portant to note that this data was taken off-resonance from the 
optical cavity (A/k:^ 1) so that the mechanical motion would 


not dominate the response; however, a trace of the mechanical 
response is still visible in the noise of the curve near the me¬ 
chanical resonance frequency (vertical grey dashed line). The 
BHD phase 0/^ was set to be sensitive to phase fluctuations of 
the reflected light, as in the experiment. Since the intensity 
modulator changes amplitude, at this homodyne phase the de¬ 
tected amplitude response is attenuated and should be consid¬ 
ered in arbitrary units. Nonetheless, in the frequency range of 
interest (0.8 -1.2 MHz) the amplitude response is qualitatively 
similar to the response of the mFALCllO. Of notable differ¬ 
ence, however, is that the overall delay time is now compara¬ 
ble to the mechanical period (see horizontal grey dashed line 
in Fig.pT[)). This delay time was measured to have the follow¬ 
ing components: 200 ns from the optical detector, 400 ns from 
the band pass Altering of the feedback circuit, 180 ns from the 
phase-shifter, and the remainder is dominated by the total sig¬ 
nal path length. The non-negligible delay in the rest of the 
feedback loop requires the inclusion of the phase shifter which 
allows us to apply a feedback force that purely damps with 
little spring shift. A consequence of this nearly constant time 
delay is that the phase of the feedback signal varies linearly 
with frequency. This becomes important when considering the 
wideband feedback loop noise in the next section. 

3. Feedback noise 

This section shows the wide span excess noise (above vac¬ 
uum noise) observed in our measurement setup at high opti¬ 
cal power. It also describes how the time delay in the feed¬ 
back loop from the previous section manifests itself as excess 
noise in the optical power spectral density at high optical pow¬ 
ers. Figure [T^ shows a wide span of the noise power spec¬ 
tral density (NPSD) of our device with the feedback engaged, 
the probe laser on resonance with the cavity (A = 0), and the 
homodyne phase set to the phase quadrature (0/j = 7r/2). At 
low optical signal power (nc < 100), the noise floor is set by 
the vacuum noise of the signal arm and the Lorentzian me¬ 
chanical response dominates the signal. At higher signal pow¬ 
ers (nc = 734, green) there is ^ 3 dB additional broadband 
noise above vacuum noise which drops off with increasing fre¬ 
quency, but the mechanical response near 940 kHz still dom¬ 
inates. This broadband noise is also seen in the open loop 
measurements of Fig. [T^. A separate measurement indicates 
that this background noise is not laser intensity noise. Ad¬ 
ditionally, by measuring the same broadband noise off reso¬ 
nance from the optical cavity (A/k 1) we have eliminated 
the cavity (through thermo-refractive effects) as the source of 
this excess broadband noise. The noise is therefore likely due 
to residual optical phase fluctuations in our set-up, possibly 
due to intrinsic phase noise of the laser or acoustic noise pick 
up . Note that in principle homodyne detection is insensitive to 
frequency fluctuations of the laser since the laser is interfered 
with itself when a common LO is used, however, in practice 
the two arms of the interferometer are never exactly the same 
length and the interference is between light emitted by the laser 
at different times resulting in a conversion of frequency noise 
to intensity noise which is detected. Currently this limits our 
imprecision to about ^imp = lO”"^ quanta. 

In the experiments presented in the main text it is not, how- 
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FIG. 10: a, Magnitude (blue) and phase (green) response of the mFALCllO for the settings used in the measurements presented in the main 
text. Vertical grey dashed line indicates the mechanical frequency and the cyan dashed line indicates the amplitude response for ideal derivative 
feedback, proportional to co. b, Total time delay (red) as measured on a network analyzer. Horizontal grey dashed line indicates one mechanical 
period. 
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FIG. 11: a, Magnitude (blue) and phase (green) response of the entire feedback system including electronics and optical train, with the laser 
off resonance from the optical cavity. The grey dashed line indicates the mechanical frequency. The noisy signal centered at the mechanical 
frequency is simply a result of residual mechanical transduction and is not indicative of a rapidly changing response of the circuit, b. Total 
time delay as measured on a network analyzer (red). A single mechanical period is indicated by the horizontal grey dashed line. 


ever, the broadband excess imprecision noise which ends up 
limiting the attainable cooling, as can be seen in Fig. l^i, at 
increased laser power {n^ = 925; red curve) there is a broad 
noise peak centered near ^ 800 kHz which grows and pre¬ 
vents further cooling (or accurate determination of the mode 
occupancy for that matter). In order to investigate this further 
we measure the NPSD with the laser off-resonance (A/k 1) 
and the feedback engaged in Fig. 13 Note that because the 
data was taken off-resonance, there is very little intracavity 
power. For direct comparison then, we have taken the data 
at = 500, which is the number of photons that would be in 


the cavity for the same reflected optical power as with the laser 
directly on resonance (A = 0). The feedback loop produces ex¬ 
cess noise “humps” at approximately 800 kHz, 1.5 MHz, and 
2.2 MHz frequencies. The “humps” in Fig, correspond in 
frequency to where the phase of the entire feedback loop (see 
Fig. [TT^) is between 0® and 180^ modulo 360^, which cor¬ 
responds to the regions where the feedback loop has positive 
gain, amplifying the measured noise rather than damping it. 
This limits the bandwidth over which we can effectively damp 
and cool the mechanical motion to about 250 kHz, which for 
the present device limits the attainable cooling level. 
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FIG. 12: a, The wide-span noise spectra for a series of experimental powers with the feedback loop engaged (closed loop) and the laser tuned 
on-resonance with the optical cavity (A = 0,0/^ = 7i/2). b, The wide-span noise spectra for a series of experimental powers with the feedback 
loop disengaged (open loop) and the laser tuned on resonance with the optical cavity (A = 0,0/^ = 7i/2). The noise floor at the lowest photon 
number shown is due to vacuum noise and its slight downward slope is from the detector’s response. Note that the discrepancy between the 
absolute noise levels between the two plots is because the data were taken at different detector gain and local oscillator powers. 
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FIG. 13: Wide-span noise spectrum for the feedback loop engaged 
(closed loop) and the laser far detuned from the laser cavity (A/k >> 
1). Here we plot the spectrum for a series of different homodyne 
detector phases around 0/^ = 7t/2. Measurements were taken at = 
500 which corresponds to the number of photons that would be in the 
cavity to yield the same reflected power at the detector if the laser 
were on resonance A = 0. 































